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Light Curves and Spectra of Dust Echoes From Gamma-Ray Bursts and their 
Afterglows: Continued Evidence that GRB 970228 is Associated with a 

Supernova 

Daniel E. Reichart 1 ' 2 
ABSTRACT 

The relative brightening and reddening of the optical afterglows of GRB 970228 
and GRB 980326 around 20 - 30 days after these bursts have been attributed to 
supernovae, which are red (blueward of their spectral peak), and which peak in time 
after ~ 20(1 + z) days. However, this direct evidence for a GRB/SN connection has 
recently been challenged. It has been suggested that the late afterglows of these bursts 
can be explained by dust echoes, of which we consider two cases: (1) the scattering 
of light from the afterglow (the forward shock), the optical flash (the reverse shock), 
and/or optical light from the burst itself by a shell of dust at a radius R from the 
progenitor, and (2) the absorption and thermal re-emission of this light by this shell 
of dust. In this paper, we model and compute dust echo light curves and spectra for 
both cases. Although the late-time afterglow of GRB 980326 was not sufficiently well 
sampled to rule out a dust echo description, we find that the late-time afterglow of 
GRB 970228 cannot be explained by a dust echo. 

Subject headings: dust, extinction — gamma-rays: bursts 

1. Introduction 

Bloom et al. (1999) first identified a second component to the optical afterglow of GRB 
980326. This second component grew brighter than the extrapolation of the early afterglow 
(which faded as a power law) around 1 week after the burst, and peaked around 20 - 30 days 
after the burst before fading away. Rough spectra taken 2.4 and 28 days after the burst show 
that the second component was significantly redder than the early afterglow. Bloom et al. (1999) 
attributed this behavior to a supernova at a redshift of z ~ 1. Supernovae peak after about 
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20(1 + z) days and are significantly redder than afterglows (blueward of the supernova's spectral 
peak). Unfortunately, neither a redshift nor multiband photometry is available for this afterglow. 

Reichart (1999) identified a similar second component to the optical afterglow of GRB 970228, 
the first burst for which an afterglow was detected. Reichart (1999) and Galama et al. (2000) 
showed that the light curve of this second component is consistent with the expectation of a 
supernova at the spectroscopically-measured redshift of this burst (z = 0.695; Bloom, Djorgovski, 
& Kulkarni 2001). More convincingly, VRIJK photometry of the second component taken near 
its peak showed that its spectral flux distribution is also consistent with the expectation of a 
supernova at this redshift: (1) the spectral flux distribution is very red blueward of the spectral 
peak and very blue redward of the spectral peak, and (2) the spectral peak is at the frequency 
of the spectral peak of a supernova at z ~ 0.7 (Figure 1; Reichart 1999; Galama et al. 2000; 
Reichart, Castander, & Lamb 2000; Reichart, Lamb, & Castander 2000). This is the best evidence 
to date for the GRB/SN connection. 

The GRB/SN connection has recently been challenged by alternative explanations of these 
observations, called dust echoes (Waxman & Draine 2000; Esin & Blandford 2000; see also Lamb 
1997, 2000). If the progenitors of bursts are massive stars (see Lamb & Reichart 2000 for a 
review of the indirect evidence that bursts are associated with star- forming regions), due to the 
progenitor's wind, the density of the surrounding medium should fall off as n(r) oc r~ 2 out to 
the radius R w of the wind termination shock (Chevalier & Li 1999). At R w , the density should 
jump by a factor of « 4, and beyond R w , the density should remain approximately constant out 
to great distances (Chevalier & Li 1999). Dust in the vicinity of the burst should be sublimated 
by UV light from the optical flash (the reverse shock; Waxman & Draine 2000). Let R be the 
radius out to which dust is sublimated. If R < R w , and if the dust is sublimated isotropically, 
the dust should occupy a relatively thin shell at ~ R due to the rapidly decreasing density of the 
circumburst medium with distance. Light that is emitted by the afterglow (the forward shock) 
from a distance that is much less than R at an angle 9 relative to the line of sight, is scattered 
(Esin & Blandford 2000) or absorbed and thermally re-emitted (Waxman & Draine 2000; Esin & 
Blandford 2000) by dust at r ps R into the line of sight, and arrives at the observer, arrives with 
an observer-frame time delay of 

t= r(1 + Z) (l-cosg) (1) 
c 

relative to light that was emitted into the line of sight (6 = 0) and not scattered or absorbed. This 
dust and light geometry is depicted in the top panel of Figure 2. Waxman & Draine (2000) and 
Esin & Blandford (2000) find that this results in secondary components to the afterglow, called 
dust echoes, that may peak about 20 - 30 days after the burst, and that may be significantly 
redder than the early afterglow, i.e., that the late-time afterglows of GRB 970228 and GRB 980326 
may in fact be due to dust echoes, and not supernovae. 

There are a number of problems with this dust and light geometry. First of all, R ~ 12 pc 
(Waxman & Draine 2000); however, R w is probably only a few parsecs (Chavalier & Li 2000). 
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Consequently, the dust would probably not occupy a thin shell. 

Secondly and more importantly, UV light from the optical flash will only sublimate dust 
within Oj e t of the line of sight, where 8j e t is the initial half-opening angle of the jet. (The Lorentz 
factor 7 of the shocked material in the jet is much greater than 6~^ t when the optical flash occurs.) 
If the shock is adiabatic, 

, „ 6 . 1(1 + ,)3/s - 1/8 (_*_) 1/8 («) - 3/8 , (2) 

' VlcmV V10 52 er g y \ldyj V ' 

where n is the number density of gas atoms in the circumburst medium, E is the isotropic- 
equivalent energy of the shock, and t is observer-frame time (e.g., Sari, Piran, & Narayan 1998). 
Let n = n(l pc)(r/l pc)~ 2 . Since r « 47 2 ci/(l + z) (e.g., Sari, Piran, & Narayan 1998), 



7 » 2.2(1 + 2 



,V4 



n(l pc) 



1 cm 



-1/4 / E \ 1/4 / + \ -1/4 

' ' ' 1 . (3) 



10 52 erg/ \ldy 



The afterglow emits light into a half-opening angle Q co \ ~ (# 2 et + 7 _2 )i/ 2 . Given the rate at which 
7 decreases with time, 6 co i quickly grows significantly greater than 9j et : For example, 9 co i > 1.16 j et 
after only a few minutes if 6j e t ~ 0.1^| 



n(lpc)" 1 1 
1 cm 3 
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Light emitted at 6 > 6j e t would be scattered, absorbed, and thermally re-emitted by dust at radii 
that are substantially smaller than R. Due to the substantially higher density of dust at these 
radii, this light would dominate the dust echo. Radiative shocks yield similar results. This more 
realistic dust and light geometry is depicted in the bottom panel of Figure 2. 

In this paper we consider a dust and light geometry of which the dust and light geometry 
depicted in the bottom panel of Figure 2 is a limiting case. High-energy light from the progenitor 
may destroy dust (e.g., via electrostatic disruption; Draine & Salpeter 1979; Waxman & Draine 
2000) out to a smaller radius. Similarly, irregularities in the progenitor's wind may destroy dust 
in shocks. In either case, the dust would be destroyed isotropically. Let R instead be the radius 
out to which dust is destroyed by the progenitor, and let Rj e t be the radius out to which dust 
is destroyed by the burst, the optical flash, and/or the afterglow, by whatever dust destruction 
mechanisms. If Rj e t is even twice as great as R, and if the jet is oriented along the line of sight, 
the dust echoes should brighten substantially at a time given by Equation (|]) with r ~ R and 



3 When 7 ~ 0. ., the jet begins to expand laterally, and the afterglow fades much more rapidly (e.g., Panaitescu, 
Meszaros, & Rees 1998). Applying Equation (||) to the afterglows of GRB 990510 (e.g., Stanek et al. 1999; Harrison 
et al. 1999), GRB 991216 (Halpern et al. 2000), GRB 000301c (e.g., Rhoads & Fruchter 2000; Berger et al. 
2000), and GRB 000926 (Price et al. 2001), we find that 9 jet « 0.17[n(l pc)/l cm 3 ] 1/4 , 0.14[n(l pc)/l cm 3 ] 1 ' 4 , 
0.38[n(lpc)/lcm 3 ] 1/4 , and 0.18[ra(l pc)/l cm 3 ] 1 ' 4 , respectively. 
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9 ~ &jet, and dim substantially at a time given by Equation (||) with r ~ i? and # ~ 7r — 9j e t, since 
we assume that the jet is bipolar. If 9j e t -C 1, then the dust echoes "turn on" at 

eo „»6.0( 1 + ,) (^) (|f) 2 dy, (5) 

and "turn off" at 

t o// « 6.5(1 + z) (— J yr - t on , (6) 

though they may fade earlier than t ff due to other effects (see §2 and §3). The jet may be 
oriented as much as ~ 9 jet from the line of sight (any more and we would not have seen the burst). 
Consequently, the dust echoes may "turn on" half of the way as early as t ~ and all of the way 
as late as t ~ 4i on . The "turn off" time would be similarly affected. The dust and light geometry 
depicted in the bottom panel of Figure 2 is recovered in the limit that R — > 0. 

Another dust and light geometry that we consider in this paper is that light from the optical 
flash, as well as optical light that may be associated with the burst itself, may result in dust 
echoes. Consequently, dust echoes may be a way to measure the brightness of the optical flash 
and/or the optical brightness of the burst. 

In this paper, we model and compute dust echo light curves and spectra for both Rj e t > a few 
times R, and Rj e t < R (in which case dust is destroyed only isotropically) . Any other case would 
be fine tuning. We also model and compute dust echoes resulting from a pulse of light (i.e., the 
optical flash or optical light associated with the burst), and an afterglow light curve. Finally, we 
model and compute dust echo light curves and spectra for both scattered light (§2), and absorbed 
and thermally re-emitted light (§3). We find that both types of dust echoes have light curves 
and spectra that can be distinguished from each other, and from the light curves and spectra of 
supernovae. In §4, we re-examine the late-time afterglow of GRB 970228, and find that it cannot 
be explained by a dust echo. We draw conclusions in §5. 



2. Dust Echoes from Scattered Light 

2.1. Dust Echoes from a Scattered Isotropic Pulse of Light 

Before modeling the general case of the dust echo of collimated light that is emitted as a 
function of time (e.g., an afterglow), we first model the simpler case of the dust echo of an isotropic 
pulse of light (e.g., an optical flash or optical light associated with the burst) that is emitted from 
a distance that is much less than R. Consider the light that is emitted at an angle 9 relative to 
the line of sight, is scattered by dust at r into the line of sight, and arrives at the observer. The 
optical depth at observer-frame frequency v along this path is given by: 
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where n(r, 9) is the number density of dust particles at r and 9, and c*°' 1+2 ) = a v(i+ z \ + a t\i+z) ^ s 
the total (absorption plus scattering) cross section of these particles. The number density of dust 
particles at r < R w is given by (§1): 



(r < R) 

. (r < R jet and < 9 < 9 jet or vr - 9 jet < 9 < vr) 

" ! ''• ' ^ n R (r/R)- 2 (R<r<R w and /0, et < < vr - /% et ) ' 1 j 

n R (r/R)- 2 (fya < r < R w \i R < R jet < R w ) 



where (§1) 



.0 (^<i?) 
1 1 (i?jet > a few times i?) 



Assuming that R w > a few times R, we find that 



T u(l+z) 

<v(i+z)(r,9) = I Ty{1+z) 



1 + T ~ !)] < < f or rsinfl > fl) 



1 + S 



fl ^-/e jet -2cos- 1 ^ 



sin 



(| < 9 < 7T - /0j ef and rsin# < R) 



Mi+zp-^) (0 = 7rand/ = O) 

(10) 

where t j/ ( 1+;S ) = urg^R is the optical depth along any line of sight if / = 0, or any line of sight 
between 9j et < 9 < it — 9j e t if / = 1. 

The probability that light takes this path, i.e., that light that is emitted within dfl of 9 makes 
it to within dr of r without being scattered or absorbed, is scattered within dr of r into the line of 
sight, and makes it to the observer without being scattered or absorbed, is given by: 

d K(i+z)( r > e ) = e- T ^+^ r ' e) n(r, 9) dUdr. (11) 

We adopt the same differential dust scattering cross section as do Esin & Blandford (2000; White 
1979): 



aa u{l+z) _ e Kl+z)°V(l+2) 1 ~~ VviX+z) 



(12) 



d& 4vr [1 + fil (1+g) - 2/x, (1+2)/ i] 3 / 2 ' 

where \i = cos 9, ji v r\+ z \ is the mean value of /x at v{l + z) for the dust particles, and 
e v(i+z) = a tli+z)/ a °(l+z) * s ^ ne mean albedo at v(l + z) for these particles. We adopt 
fivfi+z) = e u(i+z) = 0-5 (White 1979). The solid angle across which this light is scattered is 
d£l = 2tt sin 9 d9. Consequently, 

- ^" + " " Ml q^f^-i (L) . (is) 

u{1+z)y 1 2 n R (5-4cos#) 3 / 2 \RJ V ; 

Using Equation (|l]), we rewrite Equation ( |i~3| ) as a function of r and t: 

dP sc (r t) ~ 3CT ^ l+z) " M) ("Y 1 e ~ Ml+2)M) d (L\ dt ( U ) 
dP »(i+*)^ t )- 2 R(l + z) n R \RJ (5-4cos0)3/2 d UJ ' (U) 
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where 6 is a function of r and t given by Equation ([!]). 



Let F v dt be the fluence of the pulse of light at v. Ignoring light that multiply scatters into 
the line of sight, the light curve of the dust echo is then given by: 



,-(i-/K (1+z) (iP Ki+^( r ' t ) 



3cTj,( 1+ ^) 

2R(l + z) 



P(r') ( ^ 



(5 -4 cos 0) 3 / 2 



-d [ - ) F u dt, 



(15) 
(16) 



where 



and 



and where 



[r' < r(7r - f9 jet )} 
P(r') = { 1 [rfr - /% et ) < r' < r(/0, et )] 
[r'>r(/% e t)] 

i? [t< t(7T - f9 jet )} 

r(ir - fe jet ) [t > t(ir - fe jet )\ 

f 9(R) [t < t(7T - fd jet )} 

[TV-fe jet [t>t(ir- f9 jet )] ' 



r{6) 



ct 



6{r) = cos 1 



(l-cos0)(l + z) 
ct 



1 



and 



t{6) 



r(l +z) 
(1 - cos e)R(l + z\ 



(17) 

(18) 
(19) 

(20) 
(21) 

(22) 



If / = 0, or if / = 1 and 9 jet < min{0, (0/r)(r/R)} [Equations @ and (|l0|), respectively], then 
Equation is approximately given by: 



where 



A 
A 



\t sc 



A (7i) 2 {it 



-3/2 
-3/2 



4 " 3 (i) 

7/4 



-1/4' 



{t < t sc and t < t T ) 
(t sc <t<t T ) 

(t sc <t T <t) 

(t T < t and t T <C t sc ) 



(23) 



A = 3CTHl+z \ e- fT »(^F„dt, 
4R(1 +z) 

lR{l + z) ( r 
R 



(24) 
(25) 
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is the time when the scattering cross section begins to fall off due to scattering at greater angles, 
and 



3 R{l + z) 



r V 

RJ 



(26) 

T u(l+z) c 

is the time when the optical depth grows greater by 1 due to the greater path length at greater 
scattering angles. This interaction of the dust and light geometries is depicted in the top panel of 
Figure 3. We plot example light curves for different values of /, R(l + z), and r w (i +z ) in the top 
panels of Figures 4-7; however, consider the following limiting cases: 

f = (Figures 4 and 5, top panels): If t„( 1+z ) <J 12, then t sc <; t T . Consequently, F^{t) oc t° 
from t = until t ~ i sc , after which the light curve slowly rolls over to F^{t) oc i~ 3 / 2 . If 
t v [x-\- z ) ^ 1-5) then this continues until t ~ i //, after which F^{t) oc i~ 2 . If 1.5 <J t u /i +z \ <^ 12, 
then this continues until t ~ t T , after which the light curve slowly rolls over to F^(t) oc i~ 7 / 4 until 
after t = t ff. After this time, F^{t) oc i~ 2 . If t„( 1+2 ) » 12, then t T is initially less than t sc . 
Consequently, F^{t) oc i° from t « until i « i T , after which F^{t) oc i -1 until after t = t ff. 
After this time, F^{t) oc t" 2 . 

f = 1 (Figures 6 and 7, top panels): The light curve has the same dependence on time as in 
the / = case except that it does not "turn on" until ~ t on . Also, since the dust is destroyed 
along the line of sight, the pulse of light is a factor of e Tv( - 1+z '> brighter than in the / = case. 
Consequently, the light curve, which is normalized to the fluence of the pulse of light, is a factor 
of e _r,/(1 + z) fainter than in the / = case. 

Let T v n +Z \ oc [u{l + z)] Q "( 1 + z ) . The value of a u n +z \ is canonically assumed to be 1; however, 
the true value is closer to 1.6 at NIR values of v{l + z), and levels off to ss 1 in the UV (Reichart 
2000). The spectral index of the dust echo at v as a function of t is then given by: 



dlogF u E (t) 
dlog v 

a u + QV(i+z) 



(27) 



1 ~ 217 ~ 7*1/(1+*) 



10 



±(ii 



-5/2 



- 3 (^J 



f T v(l+z) 



1 - 



i gz3jferr(i) 

10 -3(^)- 1/2 (-) 



1/2 t 

5/2 ^ t ^ -5/4 

T74 — f T u(l+z) 



\ ~ fTv(l +Z ) 



(t < t sc and t < t T ) 

{t ac <t<t T ) 

{t sc <t T <t) 

{t T < t and t T <C t sc ) 



,(28) 



where a v is the spectral index of the pulse of light. This approximation holds if / = 0, or if / = 1 
and 9 j e t <C min{#, {6 /r){r / R)}. We plot how the spectral index of the dust echo differs from the 
spectral index of the pulse of light as a function of time for different values of /, R{1 + z), and 
T u(i+z) hi the bottom panels of Figures 4-7; however, consider the following limiting cases: 

f = (Figures 4 and 5, bottom panels): Since blue light scatters with a higher probability 
than red light, and since both the dust echo and the pulse of light are extinguished equally at 
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t = 0, the dust echo is initially bluer than the pulse of light, independent of the value of r u (i +z y. 
a v(t) ~ a v ~ ~ 1- If T u(i+z) & 12, then t sc <^ t T . Consequently, the dust echo reddens 

to a„(t) — a u a v n +z \(l — T v n +z \/24) at t « t sc . If t u (i +z ) ^ 1.5, then it grows redder to 
a u(t) ~ a v ~ a i^(i+^)(l ~~ t j/(i+^)/8) at t rj t // , after which it grows bluer and returns to its 
original color. If 1.5 £ r u n +z -\ <^ 12, then it grows redder to 0.5a u (x +z ) <^ o^{t) — a v ^ 0.7a„n +2 ) 
at i w t T , after which it grows bluer to a color of a„(t) — a v r* 0.8a v n +z \ until after t = t ff- After 
this time, it grows bluer still and returns to its original color. If t v (i+z) ^ 12, then t T is initially 
less than t sc . Consequently, the dust echo reddens to a^(t) — a u « 0.5a u n +z \ at t « t T , after which 
its color does not change until after i = t ff- After this time, its color grows bluer and returns to 
its original color. 

f = 1 (Figures 6 and 7, bottom panels): The spectral index of the dust echo has the same 
dependence on time as in the / = case, but is redder than the pulse of light by an additional 
T v(i+z)'- since the dust is destroyed along the line of sight, the pulse of light is not similarly 
reddened. 

In any case, if T v t\+ X \ ^ 12, then the spectral index of the dust echo does not change in time 
by more than Q.^a v (x +Z ) ~ 0.5, and if t v (\ +z ) <C 1, then it does not change at all. If t v (\ +z ) ^ 12, 
then the dust echo is probably too extinguished to be observed; consequently, we do not explore 
this case in detail in this paper. 



2.2. Dust Echoes from a Scattered Collimated Pulse of Light 



We now model the more general case of the dust echo of a collimated pulse of light that 
is emitted from a distance that is much less than R. Let the pulse of light be emitted between 
< 9 < 8 co i and ir — 9 co i < 9 < it (we assume bipolar emission), where 9 co i > 9j et . The only 
difference from the isotropic case is that dP^ 1+z ^(r, 9) = if 9 co i < 9 < ir — 9 co i. Consequently, the 
light curve is again given by Equation fll6|), but with 



P(r') 



[r> < r(ir - f9 jet )] 
[r(ir - fG je t) < r' < r(ir - 9 cot )} 
[r(7r - 9 col ) <r'< r{9 col )] 
[r{9 col ) <r'< r(f9 jet )] 
[r' > r(f9 jet )] 

[t < t(0«rf)] 

Wed) <t<t(n- 9 col )] 
[t(n - 9 col ) <t<t{ir- f9 jet ) 



r o 
i 



1 

x 

R 

r(9 col ) 
R 

r(7T - fe jet ) [t > t(lT - f9 jet ) 



(29) 



(30) 
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and 

[t < t(6 col )} 

[tiPcot) < t < t(n - 9 col )] 

[*(7T - &coi) <t<t(ir — fOj et )] ■ 1 > 

[t > t(7T - fO jet )] 

This interaction of the dust and light geometries is depicted in the bottom panel of Figure 3. 

We plot example light curves and color histories for different values of 9 co i and T v ri+ Z ) i n 
Figures 8 and 9, but basically, they are the same as in the isotropic case until t « t{6 co {). After 
this time, the light curve falls off as F^{t) oc t~ 2 if t < t T , as F^{t) oc i -7 / 4 if t > t T and t sc < t T 
[ T v(i+z) & I2r/R], and as F^(t) oc t" 1 if t > t T and t T < t sc [T u n +z -\ S> 12r/J2], because the 
observed light no longer includes light that scattered close to the evacuated region where the 
density (and the optical depth) is highest. A direct consequence of this is that the dust echo 
grows bluer. The light curve and color history partially recover their isotropic characteristics after 
t ~ t(ir — 9 coi), when back-scattered light from the far side of the evacuated region can first reach 
the observer; however, the light curve falls off rapidly after t ~ toff- 



2.3. Dust Echoes from Scattered Afterglow Light 



We now model the even more general case of the dust echo of collimated light that is emitted 
as a function of time (i.e., an afterglow) from a distance that is much less than R. Since the 
afterglow emits light continuously, at any time t after the burst, we receive scattered light that 
was emitted by the afterglow at all prior observer-frame times: < t' < t. Consequently, the dust 
echo light curve is given by integrating Equation (|i~6|) over these times: 



'\~ 3 p -[ T v (i+z)(r' ,9)-(l-f)T v(1+z) ] , I 

-d[-)F v {t')dt\ (32) 



(5 -4 cos 9fl 2 



R 



where P(r'), r, and 9 are given by Equations (|29|), (p0|), and (|3l|) , respectively, but with 

c(t - f ) 



r{9) 



9{r) = cos 



(l-cos6»)(l + z) 
c(t - t') 



-i 



and 



t(9) =t' + 



r(l + z) 
(1- cos 0)12(1 + z 



(33) 
(34) 

(35) 



F u (t') is the light curve of the afterglow, and 9 co \ is now a function of time that we approximate 
by: 

roi = min((0? t + 7 - 2 ) 1/2 ,?)- (36) 
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Let tj e t be the time when 7 = - et . Then from Equation (|3|), 
and 

^ = 5.5(1 + ^^(1 pc) 



1 cur 




(37) 



(38) 



The value of 7 oc (t') after tj e t (e.g., Sari, Piran, & Halpern 1999); this is true in either a 
n(r) oc r° or a n(r) oc r~ 2 circumburst medium. 

Finally, we model the light curve of the afterglow. In this paper, we simply take: 

(F U E ^ (f<t„) 
Fftf) = Fu (&) ^ <f< t jet ) ; (39) 

which is typical of observed afterglows after to, where to may be as little as minutes or as much as 
hours after the burst (e.g., Sari, Piran, & Narayan 1998). Before to, the light curve can depend 
on time in a variety of ways (e.g., Sari, Piran, & Narayan 1998); however, we do not model it any 
more carefully than a constant. Instead, we show that the value of to has little effect on the light 
curve and color history of the dust echo (see Figure 13). 

We plot example light curves and color histories for different values of 0j e t, tj et , and to in 
Figures 10 - 13. The results are similar to the case of a collimated pulse of light with collimation 
half-angle given by the value of 9 co \ when most of the afterglow light is emitted; i.e., by the value 
of Ojet, since most of the afterglow light is emitted between to < f < tj e t, and 9 co i ~ 9j e t during 
these times. 

In Figure 10, we take / = 1, R(l + z) = 0.1 pc, t v{1+z) = 1, t = 1 hr, 6 jet = 5°, 10°, 
20°, and 40°, and tj e t = 0.35, 1.4, 5.6, and 22 dy, respectively. These values of tj e t are given by 

g2 _ 1n 52 



Equation (|38| ) with z = 1, n(l pc) = 1 cm , E6j et = 10 erg, and these values of 9j e t- Clearly, 
for these values of the parameters, the dust echo of the afterglow is substantially fainter than the 
afterglow itself. Changing the value of t u ^ +z ^ only makes the dust echo fainter relative to the 
afterglow: decreasing r^( 1+2 ) decreases the amount of light scattered into the line of sight, and 
increasing r v n +z \ increases how much the dust echo is extinguished, but does not increase how 
much the afterglow is extinguished since / = 1. The dust echo can be made brighter than the 
afterglow if / = and T v n+ Z \ is substantially increased, but in this scenario, both the afterglow 
and the dust echo are substantially extinguished and consequently difficult to detect. This leaves 
only three scenarios in which a dust echo can be brighter than an afterglow and detectable: (1) 
the value of R(l + z) is substantially greater than 0.1 pc, in which case the dust echo peaks at a 
substantially later time (but with a substantially lower peak brightness), after the light curve of 



-li- 



the afterglow has rolled over and "gotten out of the way" of the dust echo (Figure 11); (2) the 
value of (1 + z)[n(l pc)] -1 E6j et is substantially less than the value assumed above, in which case 
the light curve of the afterglow rolls over at a substantially earlier time, again "getting out of the 
way" of the dust echo (Figure 12); and (3) the dust echo is dominated by scattered light not from 
the afterglow, but rather, from the optical flash and/or optical light from the burst itself. This 
last scenario can be modeled simply as a collimated pulse of light (§2.2) with 9 co i ^ 9j e t- We show 
that the choice of to has little effect on the light curve and color history of the dust echo of the 
afterglow in Figure 13. 

One aspect of the dust and light geometry that we do not model in this paper is that the 
afterglow is not emitted from the center of the evacuated region, but rather from the head of 
the jet, which is moving away from the center (toward the observer) at relativistic speeds. Since 
r ~ 4j 2 ct/(l + z) (§1), the head of the jet reaches the radius r at which most of the scattering is 
occurring at time 

Our model holds as long as most of the light is emitted before some reasonable fraction of t r . 
For afterglow light curves like the one described in Equation (^), most of the afterglow light is 
emitted before tj e t, by a factor of a few. Consequently, our model only holds as long as t r ^> tjet, or 



r ;> 0.08 



n(l pc) 
1 cm 3 



E9 



2 

jet 



10 52 erg 



pc. (41) 



However, if the dust echo is dominated by scattered light from the optical flash and/or optical 
light from the burst, as opposed to light from the afterglow, then the light is obviously emitted 
long before t r , and consequently, the model holds for substantially smaller values of r as well. 



(42) 



3. Dust Echoes from Absorbed and Thermally Re-Emitted Light 

The dust echo from absorbed and thermally re-emitted light is easier to model than the dust 
echo from scattered light: unlike Equation (ll2j), the differential cross section for dust absorption 
has no dependence on 9: 

da l\i+z) = i 1 -Mi+^Kq+z) 
dVt ~ 4vr 
Consequently, the light curve of the dust echo of the afterglow is given by 

^^l^iTM^ VlM,+i,, ' w (43) 

where FjFt(t) is the bolometric light curve of the dust echo, and Fuy{t') is the light curve of the 
afterglow, integrated over UV values of v{l + z). If the re-radiating dust particles had a common 
temperature T, the spectrum of the dust echo would be thermal. However, the temperature of 
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the dust particles decreases with radius as T oc r -1 / 2 (Waxman & Draine 2000); consequently, 
the spectrum should fall off more slowly than a blackbody redward of the spectral peak. Since 
dust is rapidly sublimated at temperatures T ^> 2300 K, the spectral peak should occur only at 
frequencies v <^ 1.35 x 10 14 (1 + Hz. For a burst at z = 0, this is in the K band. 

We plot example light curves for different values of 8j e t in Figure 14. However, consider the 
limiting cases for an isotropic pulse of light. The thermal dust echo has the same limiting cases as 
the scattered dust echo (e.g., Equation |23|), except that t sc = oo. Consequently, if t T > t(6 co i), then 
Fjn(t) oc t° until t ~ t(6 co i), after which Fjn(t) oc t~ 2 until t ~ t T . After this time, Fm(t) cx t . 
If t T < t(9coi), then Fm(t) oc t° until t m t T , after which Fjjt(t) oc t . 



4. Can the Late-Time Afterglow of GRB 970228 Be Explained by a Dust Echo? 

No. First of all, the late-time afterglow of GRB 970228 cannot be explained by a thermal dust 
echo, because its spectral flux distribution (Figure 1) peaks between the V and B bands in the 
source frame, not in or beyond the K band (§3). Secondly, the late-time afterglow of GRB 970228 
cannot be explained by a dust echo due to scattered afterglow light, because its light curve (e.g., 
Reichart 1999) does not roll over at sufficiently early times for the dust echo to be able to outshine 
the afterglow when the second component to the light curve peaks, about 20 - 30 days after the 
burst (§2.3). This leaves only the possibility of a dust echo due to scattered light from the optical 
flash and/or optical light from the burst itself. To investigate this possibility, we have fitted the 
spectral flux distribution of Figure 1 with two power laws, one to the reddest three points and 
one to the bluest three points. We find that the spectral index redward of and including the 
source- frame V band is af v = 2.29^ ; 5g , and that the spectral index blueward of and including 
the source-frame B band is af B = — 3.09^q;55. The corresponding optical depths as a function of 
the spectral index of the pre-scattered, unextinguished light are given by Equation (|28| ) with / set 
to 1. Taking a r y = 1.6 and a^s = 1 (§2-1), we find that for an intrinsic spectral index of -1/2, 
r r y = — 0.99^Qg3 and t^b = 3.34+g'go, and for an intrinsic spectral index of 1/3, r r y m — 0.47j^gjj 
and TbB = 4.17^q'qq. Intrinsic spectral indexes of -1/2 and 1/3 are typical of what one expects for 
the optical flash blueward and redward of the synchrotron peak, respectively, which one expects 
to be at optical/NIR wavelengths during the optical flash (e.g., Sari & Piran 1999). In the case 
of optical light from the burst itself, one expects a spectral index of 1/3 across the optical and 
NIR bands (e.g., Sari & Piran 1999). Consequently, in either interpretation, the value of r r y is 
borderline consistent with being physical (i.e., non-negative), but the rapid transition between Ar 
or Ay rs to Ab or Ajj sa 3.6 ± 0.7 or 4.5 ± 0.7 is unphysical. Consequently, we find that the 
late-time afterglow of GRB 970228 cannot be explained by a dust echo. 

Unfortunately, the late-time afterglow of GRB 980326 was not sufficiently well sampled to 
distinguish between a supernova and a dust echo explanation 
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5. Discussion and Conclusions 

We have modeled and computed dust echo light curves and spectra for dust and light 
geometries that are more realistic than had previously been considered. In particular, bursts are 
very likely to be collimated (on both theory and energetics grounds), and dust may be destroyed 
by the progenitor, as well as the burst, optical flash, and/or afterglow. Also, the optical flash 
and/or optical light from the burst itself may result in dust echoes. Consequently, dust echoes 
may be a way to measure the brightness of the optical flash and/or the optical brightness of the 
burst. 

We find that dust echo spectra can be used to distinguish different types of dust echoes from 
each other, and from supernovae. For example, thermal dust echoes peak in the source-frame K 
band, while supernovae peak in the source-frame V or B bands. In the absence of extinction (e.g., 
in the source-frame NIR), dust echoes due to scattered light are bluer than the pre-scattered light 
by about 1 - 1.6 in the spectral index, and consequently, can peak in the source-frame V or B 
bands, or other nearby bands, depending on the amount of extinction. However, dust echoes due 
to scattered light do not peak as sharply as do thermal dust echoes or supernovae. It is for these 
reasons that we find that the late-time afterglow of GRB 970228 cannot be explained by a dust 
echo. 

Furthermore, dust echo light curves can often be used to distinguish dust echoes from 
supernovae. While dust echoes can peak about 20(l+z) days after a burst, they also can peak 
substantially earlier and substantially later. However, supernovae can peak only at this time. 

Finally, dust echoes can be used to probe the circumburst environment, and the effects of the 
progenitor, burst, optical flash, and/or afterglow on this environment. For example, if the dust 
echo light curve brightens, levels off for a long time, and then fades, then / = and dust had 
not (yet) been destroyed along the line of sight. If the dust echo light curve brightens and fades 
without leveling off, then / = 1 and dust had been destroyed along the line of sight. The time at 
which the dust echo light curve begins to fade gives you (most likely) a degenerate combination of 
R and 9j e t- The color of the dust echo, particularly near peak, gives you t u ( 1+z ), and consequently 
a degenerate combination of R and n(l pc). The time at which the afterglow light curve rolls 
over due to lateral expansion of the jet gives you a degenerate combination of 9j e t and n(l pc). 
Consequently, R, 6j e t, and n(l pc) can be separated and solved for. Furthermore, the brightness 
of the dust echo, particularly near peak, gives you a degenerate combination of R, t„m +z ), and 
the fluence of the pre-scattered, unextinguished light; consequently, this fluence can be solved for 
as well. Also, by comparing t v i 1+z \ to the same quantity measured from the afterglow, and to 
absorption measured in X rays, one can further probe dust destruction along the line of sight, and 
the metals-to-dust ratio off the line of sight where dust is not destroyed, respectively. 

Support for this work was provided by NASA through the Hubble Fellowship grant #HST- 
SF-01133.01-A from the Space Telescope Science Institute, which is operated by the Association 
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Fig. 1. — The K- through V-band spectral flux distribution of the late-time afterglow of GRB 970228 
after subtracting out the contribution of the host galaxy and correcting for Galactic extinction (filled 
circles and upper limit), and the I- through U-band spectral flux distribution of SN 1998bw after 
transforming to the redshift of GRB 970228, z = 0.695 (Bloom, Djorgovski, & Kulkarni 2001), and 
correcting for Galactic extinction (unfilled circles; Reichart, Castander, &; Lamb 2000; Reichart, 
Lamb, & Castander 2000). The K-band upper limits are 2 and 3 a. We have fitted power laws 
to the K-, J-, and I-band observations, and the I-, R-, and V-band observations of the late-time 
afterglow of GRB 970228 (see §4). The solid lines are the best fits, and the dotted lines are the 1 a 
uncertainties in these fits. Clearly, the spectral flux distribution of the late-time afterglow of GRB 
970228 is very blue redward of the peak, and very red blueward of the peak. 

Fig. 2. — Dust and light geometries. Top panel: Isotropic dust destruction and an isotropic pulse 
of light. Dust is destroyed isotropically about the progenitor, and interior to the thick solid circle, 
the radius of which is R < R w , where R w is the radius of the wind termination shock. Exterior 
to R, n(r) oc r~ 2 out to R w . The parabola marks the surface on which light must either scatter 
or be absorbed and thermally re-emitted to arrive at an observer to the far right at a constant 
time t after the burst, where t is given by Equation (Q). Light cannot scatter or be absorbed and 
thermally re-emitted on the dotted portion of this surface of revolution because dust has been 
destroyed here. Bottom panel: Collimated dust destruction and a collimated pulse of light. Dust 
is destroyed within the bipolar cone (thick solid curve) out to a radius R > R w . The radius R w 
is marked with a thin dashed circle. Interior to R w , n(r) oc r~ 2 ; exterior to R w , n(r) oc r . The 
shaded area marks the region not illuminated by the collimated pulse of light. The parabola again 
marks a "surface of constant time" . Light cannot scatter or be absorbed and thermally re-emitted 
on the dotted portions of this surface because either dust has been destroyed here, or light does not 
reach here. In this dust and light geometry, light clearly scatters and is absorbed and thermally 
re-emitted at radii r<E. 

Fig. 3. — In both panels, the dust geometry is the same: dust is destroyed interior to a sphere of 
radius R < R w by the progenitor, and interior to a bipolar cone by the burst, optical flash, and/or 
afterglow (thick solid curve). Exterior to the thick solid curve, n(r) oc r~ 2 out to R w . In the top 
panel, the pulse of light is isotropic; in the bottom panel, the pulse of light is collimated. The 
parabolas mark surfaces of constant time that are separated by factors of two in time. Clearly, 
these dust and light geometries would result in dust echoes with similar early-time and late-time 
behaviors, but different intermediate-time behaviors. 

Fig. 4. — Light curves (top panel) and color histories (bottom panel) of scattering dust echoes of 
an isotropic pulse of light, with / = (dust is destroyed isotropically by the progenitor, but not 
additionally by the burst, optical flash, and/or afterglow along the line of sight), R(l + z) = 0.1, 
0.3, 1, 3, and 10 pc, and T v ri +e \ = 0.1. Since / = 0, the dust echo is bluer than the pulse of light. 
Since T v ri +Z \ <C 1, the color of the dust echo is relatively constant with time. 

Fig. 5. — Light curves (top panel) and color histories (bottom panel) of scattering dust echoes of 
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an isotropic pulse of light, with / = 0, R(l + z) = 0.1 pc, and T v ri +z -\ = 0.1, 0.3, 1, 3, 10. Since 
/ = 0, the dust echo is bluer than the pulse of light. However, when T u r i+Z \ p> 1, the dust echo does 
grow somewhat redder when the dust echo begins to fade away. This is due to greater extinction, 
due to the greater path length that the light must travel through the dust at the larger scattering 
angles. These larger scattering angles come into play only at late times, and are in fact responsible 
for the onset of the fading. 

Fig. 6. — Light curves (top panel) and color histories (bottom panel) of scattering dust echoes of an 
isotropic pulse of light, with / = 1 (dust is destroyed both isotropically out to R by the progenitor, 
and within 9j e t of the line of sight out to r > a few times R by the burst, optical flash, and/or 
afterglow), R(l + z) = 0.1, 0.3, 1, 3, and 10 pc, t v ii +z \ = 0.1, and 9j e t = 0.1. The light curves 
"turn on" once 9 > 6j et . Since tw 1+z ) <C 1, the color of the dust echo is relatively constant with 
time. 

Fig. 7. — Light curves (top panel) and color histories (bottom panel) of scattering dust echoes of 
an isotropic pulse of light, with / = 1, R(l + z) = 0.1 pc, r^( 1+2 ) = 0.1, 0.3, 1, 3, 10, and 9j et = 0.1. 
When Tj,(i +;2 ) > 1, the dust echo is generally fainter and redder relative to the pulse of light, because 
the dust echo is extinguished, but the pulse of light, which travels along the dust-evacuated line of 
sight, is not. However, regardless of the value of T u ri +Z \, the dust echo is bright initially, when 9 
is only slightly greater than 9j e t, and again when 9 is only slightly less than ir — 9j e t, because at 
these angles, scattered light travels through a minimal amount of dust before entering the dust- 
evacuated cone (of half-angle 9j e t ) along the line of sight. For the same reason, the dust echo grows 
significantly bluer at these times. 

Fig. 8. — Light curves (top panel) and color histories (bottom panel) of scattering dust echoes of 
a collimated pulse of light, with / = 1, R(l + z) = 0.1 pc, r v (i +z ^ = 0.1, 9j e t = 0.1 (5.7°), and 
9 col = 10°, 30°, 50°, 70°, and 90°. The 9 co i = 90° light curve and color history is the same as in 
Figure 7. The other light curves and color histories are the same as this one, except that the light 
that scattered at 9 > 9 co i no longer contributes. 

Fig. 9. — Same as Figure 8, but with t„( 1+z ) = 3. 

Fig. 10. — Light curves (top panel) and color histories (bottom panel) of scattering dust echoes 
(solid lines) of afterglow light (dotted lines), with / = 1, i?(l + z) = 0.1 pc, T v fx +Z \ = 1, to = 1 hr, 
9 jet = 5°, 10°, 20°, and 40°, and [(l+z)/2][ra(l pc)/l cm 3 ]- 1 (E9f et /W 52 erg) = 1. The corresponding 
values of tj e t are 0.35, 1.4, 5.6, and 22 dy, respectively (Equation |38|). Clearly, the dust echoes of 
these afterglows would be lost in the brightness of the afterglow. 

Fig. 11. — Same as Figure 10, but with R(l + z) = 1 pc. The dust echoes are brighter than the 
afterglows at times substantially later than tjet- 

Fig. 12.— Same as Figure 10, but with [(1 + z)/2][n(l pc)/l cm 3 ]- 1 (E^ et /10 52 erg) = 0.1. The 
corresponding values of tj e t are 0.035, 0.14, 0.56, and 2.2 dy, respectively. The dust echoes are 
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again brighter than the afterglows at times substantially later than i jet . 

Fig. 13. — Light curves (top panel) and color histories (bottom panel) of scattering dust echoes 
(solid lines) of afterglow light (dotted lines), with / = 1, R(l + z) = 1 pc, t v (i +z \ = 3, to = 0.3, 
1, and 3 hr, 9 jet = 10°, and [(1 + z)/2][n(l pc)/l cm^-^Eflj^/lO 52 erg) = 0.1. The corresponding 
value of tj et is 0.14 dy. Clearly, the choice of to has little affect on the light curve and color history 
of the dust echo. 

Fig. 14. — Light curves of thermal dust echoes (solid lines) of afterglow light (dotted lines), 
with / = 1, 12(1 + z) = 1, r„(i +J? ) = 3, t = 1 hr, 6 jet = 5°, 10°, 20°, and 40°, and 
[(1 + z)/2][ra(lpc)/lcm 3 ] -1 (E0 2 et /lO 52 erg) = 0.1. The corresponding values of tj e t are 0.035, 
0.14, 0.56, and 2.2 dy, respectively. Thermal dust echoes are relatively brighter at late times than 
scattering dust echoes, because large-angle changes in the direction of the light are just as probable 
as small-angle changes, but this is not the case with scattering. Large-angle changes are especially 
important at late times. 
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